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Abstract 

The two-particle short-range correlation functions in rapidity, az¬ 
imuthal angle and transverse momentum, following from the decay of 
statistical clusters are evaluated and discussed. 


1 Introduction 

Studies of correlations between particles produced in high-energy collisions is 
a well-known method to investigate the dynamics of the production process 
[1], They are conveniently divided into ’’short-range” when the momenta of 
the studied particles are close to each other and ’’long-range”, extending over 
large distances in momentum space. 

Already in early seventies, studies of the short-range correlations in rapid¬ 
ity led to the discovery that particle production proceeds through production 

*Address: Lojasiewicza 11, 30-348 Krakow, Poland; e-mail: bialas@th.if.uj.edu.pl 
^Address: al. Mickiewicza 30, 30-059 Krakow, Poland; e-mail: bzdak@fis.agh.edu.pl 


1 



of ’’clusters” [2], This was later confirmed by more detailed studies in other 
variables [ 1 ], although the very nature of these clusters is unclear even now. 
The problem is of importance because it touches the mechanism of hadroni- 
sation, i.e., transition from the parton system, created at the early stage of 
the collision, into the produced hadrons. This apparently non-perturbative 
transition, essential to derive the structure of the produced parton system 
from the observed hadrons, cannot be easily treated by theory. Thus a phe¬ 
nomenological analysis is needed. 

An interesting approach to this problem was formulated in the statistical 
cluster model [3, 4] which assumes that the transition from the early state 
of the process of particle production, dominated by parton interactions, pro¬ 
ceeds through an intermediate stage of clusters emitting (isotropically) the 
hnal hadrons according to the rules of statistical physics^. 

The decay distribution of such a statistical cluster at rest is taken in the 
form of the Boltzmann distribution, which for a cluster moving with the 
four-velocity becomes 

dm{p] u) ~ ( 1 ) 

where fd = 1/T is the inverse cluster temperature and p± and y are the 
transverse momentum and rapidity of the final particle. 

Although the model was originally constructed for description of the 
’’soft” processes (involving only small transverse momenta), it remains an 
interesting question to what extent it is also applicable to semi-hard and 
perhaps even hard collisions. Indeed, the patron-hadron transition being a 
soft process, happening at the very end of the parton cascade, may very well 
be universal, i.e., (quasi)independent of the mechanism of parton production. 
If this is actually the case, the statistical clusters should be visible also at 
higher transverse momenta and perhaps even in all processes of particle pro¬ 
duction at high energies. This attractive possibility was recently supported 
by the evaluation [7] of the transverse momentum spectrum of the produced 
charged hadrons. It turned out that if the distribution of the cluster trans¬ 
verse Lorentz factor 7 j_ ( 7 ^ = l+u\, where u±_ is the transverse component of 
the cluster four-velocity) follows a simple power law ~ then, surprisingly 

^ This is an attractive modification of the standard statistical model, because it allows 
to explain the observed anisotropy of the momentum spectra of particles produced in 
high-energy collisions (which is a difficulty for the statistical model), while keeping, at the 
same time, its successes in description of particle abundances [5,6]. 
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enough, the transverse momentum distribution of the emitted particles 

~ [ ^Ko{/3j±m_L)Io{/3u±p±) ( 2 ) 

dpi J 7l 

closely resembles the Tsallis formula [8] which, as is well-known [9-11], closely 
resembles the shape of the data [12-15]. See [7] for more details and [16] for 
further discussion. 

This apparent success of the concept of the statistical cluster invites one 
to study its other consequences, particularly those which may provide more 
demanding tests of the idea. Following this route, in the present paper we 
discuss the two-particle correlations and show that, indeed, they give strong 
constraints on the model and, eventually, can be even used to pin down 
possible inter-cluster correlations^. 

To determine T and k, the only free parameters of the model, we have 
htted the cluster formula (2) to the transverse momentum distribution of 
pions and kaons produced in 2.76 TeV proton-proton (p-l-p) collisions [17]. 
The £t gives k = 5 and the cluster temperature T = 140 MeV. The result 
is shown in Fig. 1, where one sees that the model reproduces the data 
with better than 20% accuracy, which is good enough for our purpose. It 
is also remarkable that the same power-law distribution (without change 
of normalization) describes well both pion and kaon distributions. We also 
checked that the model describes the charged particle spectra up to p± = 200 
GeV in p-|-p at y/s = 7 TeV [18]. 

In the next section the general formula for the two-particle correlations 
in the statistical cluster model is written down, and correlations in rapidity, 
azimuthal angle and in transverse momentum are derived. Our results are 
described in Section 3. Summary and comments are given in the last section. 

2 Two-particle correlations 

The two-particle distribution is the sum of the contribution from one clus¬ 
ter and that from two different clusters dN 2 = dN^^^ + dN!^^\ Ignoring 
correlations in cluster decay (see section 4 for further discussion) we have 

dN!^''\pi,P2) = j duW{u)dNi{pi,u)dNi{p2,u) (3) 

^To our knowledge, the first application of the statistical model for description of cluster 
decays was proposed by Hayot, Henyey and Le Bellac [2]. 
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Figure 1: The single particle distributions for pious and kaons as measured 
by ALICE in p+p collisions at y/s = 2.76 TeV, compared to the statistical 
cluster model with T = 140 MeV and k = 5. The ratios data/model are 
shown at the bottom. 

and 

= J dui J du2W{ui,U2)dNi{pi,ui)dNi{p2]U2), (4) 

where W{u) is the distribution of (four)velocity of a cluster and W{ui^U 2 ) is 
the corresponding distribution of two clusters^. For the two particle correla¬ 
tion function 

C{pi,P2)d‘^Pi±dyi(fp2±dy2 = dN2{pi,P2) - dNi{pi)dNi{p2), (5) 

^fW{u)du = (N); f duidu 2 W(U 1 U 2 ) = {N{N — 1)) where N denotes number of 
clusters. 
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with dNi{p) = J duW{u)dNi{p,u), we thus have 

C{pi,P2)d‘^Pi± dyi d?p2± dy2 = 

= dN!^''\pi,P2) + / dui / du2Cu{ui,U2)dNi{pi]Ui)dNi{p2]U2). ( 6 ) 


where 


a(ui, U 2 ) = U 2 ) - W{u^)W{u2) (7) 

is the two-cluster correlation function. 

If clusters are independent, i.e., C„(mi,M 2 ) = 0, we have 

C{pi, p2)d‘^Pi_i_dyid'^P2±dy2 = dN^^''\pi,p2). (8) 

Consider a cluster^ at rapidity Y moving in the transverse direction with the 
velocity We have 

Uq = 7 j_ cosh Y; = 'y± sinh Y; u± = 'yv±] = tanh Y. (9) 

and the formula (1) becomes 

dNi{p u) = e-/3T'xmxcosh(y-y)+/3p^uxcos(</)„-(/)) 

where 0^ and 0 are the azimuthal angles of the cluster and of the produced 
particle, respectively. 

Following [7] we take 

W{u)du ~ (11) 

where G{Y) is the distribution of clusters in rapidity.^ 


2.1 Correlation in rapidity and azimuthal angle 

We start with the correlations in rapidity and azimuthal angle. Using the 
formulae of the previous section we have 

dN^^"^ = dyid‘^Pi±dy2d‘^P2± J'y±'^d'y± J d(j)udYG{Y) 

g-/37x[mix cosh(3/i-y)+m2_L cosh(y2-y)]g/3«x[pix cos(fli„-0i)-|-p2X cos{(pu-<t>2)] ^]_2) 

"‘Henceforth we assume that clusters are uncorrelated. 

®We note that our results do not depend on the specific shape of G{Y). 
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To obtain the distributions of yi — y 2 = Ay and 0i — 02 = A(j) at fixed pij_ 
and p 2 ± we integrate over = yi + ?/ 2 , 0+ = 0i + 02 and Y. The result 

is 

CiAy, A0) ~ f [07±^™(A|/)] /o[0naI^p(A0)], (13) 

J d_L 

where Jq and are the modihed Bessel functions of the hrst and second 
kind, respectively, and 

Dm{,Ay) = sjm\^ + m|_L + 2mi±m2± cosh(A?/), 

Dp{A(j)) = + 2pi^p2± cos(A0). (14) 

Correlations in A0 at hxed pij_ and p 2 ± can be obtained by integrating 
independently yi and y 2 with the result® 

C(A0) ~ f ^^Ko[f3'y±mi±]Ko[/3'y±m2±]Io [0M±-Dp(A0)]. (15) 

J 'y _L 

For fixed Ay we have, similarly, 

C{Ay) [ yy^Ko[/3-f±D^{Ay)]Io[Pu±pi±]Io[Pu±p2±\. (16) 

J 7l 


2.2 Correlation in transverse momentum 


It is also interesting to consider the distribution of moduli of transverse mo¬ 
menta [pi±,p 2 ±]- Integrating C{pi,p 2 ) over ?/i, 1 / 2 ,0i, 02 V, 0„, one obtains 


C{pi±,p2±) ~ 


d'fx 

ll 


Ko [(3'y±mi±]Ko [07±m2±]/o [(3u±pi±] Iq [(3u±p2±]. 


(17) 


3 Results 

As explained in the introduction, the two parameters of the model: the tem¬ 
perature in the cluster decay T = 140 MeV and the power k = 5 were deter¬ 
mined from the £t to the pion and kaon single-particle transverse momentum 
distributions measured by ALICE [17]. 

®We skip the factors [27r — A0] (for A0 < tt) and [A0] (for A0 > tt) since they are 
canceled when C'(A0) is divided by the distribution of mixed events. 
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Figure 2: The two-particle correlation function (16) from a statistical cluster 
for pairs of pions and kaons with various values of the transverse momenta, 
plotted vs. /S.y = yi — 1 / 2 , the rapidity separation between the two particles. 
T = 140 MeV and k = 5. C is scaled to 1 at Ay = 0. 

The correlation functions in azimuthal angle and rapidity, given by (15) 
and (16), are shown in Figs. 2, 3. In Fig. 2 the correlation function 
C{Ay;pi±,p 2 ±) (normalised to 1 at Ay = 0) is plotted vs Ay = ||/i — 1 / 2 I, 
for pairs of pions and kaons, at various values of the transverse momenta. 
One sees that C gets narrower with increasing pij_ and p 2 ± and there is also 
some mass dependence. The correlation function C{A(();pi±,p 2 ±) is plotted 
in Fig. 3. Similar features are also seen, except that the dependence on 
particle mass is more pronounced. 

Numerical calculations show that for sufficiently high transverse momenta 
(above ~ 2 GeV) and vanishing particle masses, the two-particle correlation 
functions in rapidity and in azimuthal angle, can be approximated by Gaus- 
sians with the width squared proportional to and inversely proportional 
to the product pi±p 2 ±- The proportionality factor is close to 2k. 

Recently the GMS Gollaboration published [19] extensive studies of the 
two-particle azimuthal correlation functions in p-l-Pb collisions at ^/s = 5.02 
TeV. In Fig. 4 they are compared with the predictions of the statistical 
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Figure 3: The two-particle correlation function (15) from a statistical cluster 
for pairs of pions and kaons and with various values of the transverse mo¬ 
menta, plotted vs. A(j) = 01—02, the relative azimuthal angle between the 
two particles. In this calculation T = 140 MeV and k = 5. C is scaled to 1 
at A0 = 0. 


cluster model. One sees that the data are reasonably close to the model 
predictions at transverse momenta in the region of 1 — 2 GeV. At higher 
transverse momenta the model gives correlation functions which seem some¬ 
what too narrow." 

This result agrees with the idea that the model is applicable only at 
relatively small transverse momenta. On the other hand it may perhaps also 
indicate the presence of cluster-cluster correlations (described by the second 
term in the R.H.S. of (6)) at transverse momenta above 2 GeV. Possible 
resolution of this dilemma would require more sophisticated studies and goes 
beyond the scope of this paper. 

It would be also interesting to study balance functions [20-22]. This is 
not possible at the present stage of the model since it requires additional 

^The published CMS data [19] are modified by other physical effects, e.g., flow in p-l-Pb, 
the back-to-back peak in A(f> (which is not considered in the present paper), and by the 
procedure of the background removal. 
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Figure 4: The two-particle correlation function from a statistical cluster for 
low (left) and high (right) transverse momenta compared with the CMS 
data [19] on the short-range azimuthal correlation function measured in p-l-Pb 
collisions at ^/s = 5.02 TeV. T = 140 MeV, k = 5. All functions are scaled 
to 1 at A0 = 0. 


information on the distribution of the cluster charges. 

The transverse momentum correlation (17), divided by the product of the 
two single-particle distributions®, 


c{pi±,P2±) 


C{pi±,P2±) 

Ni{pi±)Ni{p2±) 


(18) 


is shown in Figs. 5 and 6. c{pi±,p 2 ±), normalized to 1 at pi± = p 2 ±, is 
plotted in Fig. 5 vs. the ratio \pi± —p2±\/pi± (with p 2 ± < pi±) for various 
values of pi± = 1,1.5, 5 GeV.® In Fig. 6 the value of c at p_L = pi± = p 2 ± is 
plotted vs. p±. One sees a rather fast increase of c with increasing p^. 


4 Summary and comments 

In summary, we have constructed the two-particle correlation functions in¬ 
duced by the decay of a statistical cluster. The explicit formulae were given 

®It is convenient to use this definition of c(pi_l,P 2 _l) since it is proportional to the 
number of pairs divided by the number of pairs in mixed events. 

®We checked that above pij_ = 5 GeV the curves practically do not change any more. 
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Figure 5: The transverse momentum correlation function (18) from a sta¬ 
tistical cluster for pairs of pions and kaons vs. {pi± — P 2 ±)/pi±, for various 
values of pij_. Here T = 140 MeV and k = 5. c is scaled to 1 at pij_ —p 2 ± = 0. 

for correlations in rapidity, azimuthal angle and in transverse momentum. 
Using the parameters of the model determined from the £t to the single¬ 
particle distributions {k = 5, T = 140 MeV) the correlation functions were 
evaluated. Qualitative comparison with the CMS data on azimuthal correla¬ 
tions in p-|-Pb collisions at ^/s = 5.02 TeV shows that the model works well 
at transverse momenta around 1 — 2 GeV. For larger transverse momenta the 
evaluated correlation function looks somewhat too narrow, possibly indicat¬ 
ing presence of additional inter-cluster correlations. 

Several comments are in order. 

(i) It should be emphasised that the present calculation follows directly 
from the statistical cluster model and thus contains no free parameters: the 
value of the freeze-out temperature (T ~ 140 MeV) and the parameter k = 5 
were determined from the single-particle transverse momentum distribution 
of pions and kaons. It seems remarkable that spectra of both pions and kaons 
can be described simultaneously with exactly the same cluster distribution. 

(ii) As already explained in the Introduction, we hypothesize that at the 
hnal stage of the production process, the statistical clusters are formed and 
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Figure 6: The transverse momentum correlation function (18) from a statis¬ 
tical cluster for pairs of pions and kaons vs. p± = pi± = P 2 ±- T = 140 MeV 
and K = 5. 

decay into observed particles. At high transverse momenta, jet physics is 
expected to induce correlations between clusters and thus additional corre¬ 
lations between produced particles. Detailed experimental investigation of 
this region could therefore verify universality of the cluster hypothesis and 
may also give useful information on the structure of jets. 

(iii) It would be also most interesting to measure and compare the short- 
range correlation functions in p-|-p and e’''e“ collisions in order to test uni¬ 
versality of the statistical cluster picture of particle production. Also mea¬ 
surement of correlations for various pairs of particles can be very useful in 
this respect. 

(iv) In our calculations we have ignored the correlations which may appear 
in the cluster decay. The statistical clusters are rather special objects and 
their physical interpretation, and consequently the nature of their internal 
correlations, is (in our opinion) not clear. For example, it is not obvious 
if clusters, representing multi-particles states, have a well-determined mass. 
The detailed comparison of the model with data in p-|-p or e’''e“ collisions 
should shed more light on these questions. 
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(v) In the present paper we have not discussed the baryon production, as it 
is not clear if at LHC energies the statistical model can describe correctly the 
baryon multiplicities. Within the statistical cluster model one may overcome 
this difficulty, e.g., by postulating that the clusters emitting baryons are of 
different nature than those producing mesons only. This problem is under 
investigation. 
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